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Abstract. We present a study of the early build-up of dust in high redshift galaxies. The study is based on the 
analysis of 38 Damped Lya systems (DLAs) for which we derive the fraction of iron atoms in dust form, /F e . 
The sample is representative of metal-poor galaxies in the redshift range 0.6 < z < 3.4 selected on the basis of 
their absorption Hi column density (JV(Hi) > 2 x 10 20 atoms cm -2 ). We find that the dust fraction increases 
with metallicity, from /p G ~ at [Fe/H]~ —2 dex, up to /p e — 0.9 at solar metallicity; the increase is fast below 
[Fe/H] ~ —1 dex and mild at higher metallicities. We also find some evidence for an increase of /Fe with cosmic 
time; a large fraction of the systems younger than ~ 3 Gyr has /p e 0.5. These results indicate the dust-to-metal 
ratio increases in the course of chemical evolution, at variance with the hypothesis of an approximately constant 
dust-to-metal ratio, commonly adopted in models of galactic evolution. This hypothesis is consistent with local 
and high-redshift data only when the metallicity is relatively high ([Fe/H]Jj —1 dex). The results of this work 
suggest that the main mechanisms of dust formation may be rather sensitive to the level of metallicity attained 
by a galaxy in the course of its chemical evolution. A metallicity-dependent dust production by SNe II seems to 
be the most promising mechanism for explaining the rise of /F e at [Fe/H] < — 1 dex. 
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1. Introduction 

Interstellar dust plays a key role in a variety of astrophys- 
ical processes relevant to galactic evolution and affects, 
with its presence, many measurements of basic physical 
properties of stars and galaxies. For these reasons, under- 
standing the process of dust formation and evolution on 
a galactic scale is of great interest in several areas of as- 
trophysical research. In the course of galactic evolution 
the abundance of dust is expected to track that of the 
metals synthesized and ejected from the stars. Therefore, 
the dust-to-gas ratio, T>, is expected to correlate with the 
metallicity, Z or, in other words, the fraction of heavy 
elements in dust form, /d i=a T>/Z, is expected to stay 
approximately constant (e.g. Franco & Cox 1986). In the 
local universe the existence of a Z — T> correlation is sup- 
ported by empirical evidence. For instance, Issa et al. 
(1990) found a good correlation for 6 spirals of the Local 
Group and the Magellanic Clouds, estimating V from the 
visual extinction per hydrogen atom. Schmidt & Boiler 
(1993) found a similar correlation for 23 local dwarf galax- 
ies estimating V from IRAS dust mass and H I integrated- 
emission mass. From a similar analysis Lisenfeld & Ferrara 
(1998) found a good correlation for 28 dwarf irregular 
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galaxies, but not for 16 blue compact dwarfs; they at- 
tribute these differences to variations of the star-formation 
and mass-outflow rates between these different types of 
galaxies. 

The empirical Z — T> correlation found in the local uni- 
verse is, to date, the main observational input to models of 
galactic evolution which incorporate the dust component. 
A self-consistent treatment of the metal and dust content 
in the models requires taking into account the processes 
of grain formation, accretion and destruction (e.g. Dwek 
& Scalo 1980), increasing in such a way the number of 
the parameters with few, if any, observational constraints. 
The detailed treatment of all these processes can be per- 
formed using numerical models. Following this approach, 
and requiring that accretion and destruction time-scales 
evolve in a similar fashion, Dwek (1998) found that the in- 
terstellar T> I Z mass ratio is approximately constant in the 
Milky Way. More elementary, analytical models have been 
recently proposed to reproduce the dust cycle in galax- 
ies (e.g. Hirashita 1999, Edmunds 2001, Inoue 2003). An 
approximately constant T> j Z mass ratio is generally pre- 
dicted by these studies, with the exception of Inoue (2003), 
who proposes models with increasing T> / Z in the course of 
evolution. The observed linear Z — V relation in the local 
universe is interpreted by Inoue as a sequence of constant 
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galactic age. Measuring the dust-to-metal ratios in galax- 
ies caught in their early stages of evolution is fundamental 
to probe the validity of these models. This can be done by 
studying the dust content of high-redshift galaxies, a field 
of research that is receiving a special attention given its 
potential impact on our understanding of the early uni- 
verse (e.g. Silva et al. 1998, Hirashita & Ferrara 2002, 
Morgan & Edmunds 2003). 

Quasar absorption-line systems can be used to inves- 
tigate the relation between dust and metallicity at high 
rcdshift on observational grounds. This is true for the 
Damped Ly a (DLA) systems, the class of absorbers with 
the highest HI column density (AT(HI) > 10 20 3 atoms 
cm -2 ), which are believed to originate in galaxies lo- 
cated along the quasar line of sight (Wolfe et al. 1986). 
Measuring the metallicity and the dust abundance of DLA 
systems makes it possible to probe the V — Z relation 
in primeval galaxies. The metallicity is determined from 
column-density measurements, which can be rather accu- 
rate (e.g. Lu et al. 1996, Prochaska & Wolfe 1999, Molaro 
et al. 2000, Ledoux et al. 2003). Element abundances do 
not suffer from severe ionization corrections in DLA sys- 
tems (e.g. Vladilo et al. 2001), and the metallicity level 
can be estimated pretty well, at least for elements not af- 
fected by dust depletion (Pettini et al. 1999). Measuring 
the dust abundance, however, is less straightforward. 

In principle, the amount of dust in DLA systems can be 
estimated from the reddening of the background quasar. 
This technique, however, cannot be applied to individ- 
ual systems given the uncertain, variable continuum of 
the quasar. As a consequence, quasar extinction has only 
been estimated statistically, by comparing samples with 
and without foreground DLAs (Pei et al. 1991). The dust- 
to-gas ratio in DLA systems estimated in this way is be- 
tween 5% and 20% of the Galactic value, in broad agree- 
ment with the low level of metallicity of DLA systems. 
However, this type of analysis does not give information 
on the evolution of the dust abundance with the metallic- 
ity. Also the dust emission properties cannot be employed 
to derive the dust abundance in DLA systems. In fact, 
attempts to identify the emission of intervening galaxies 
in the quasar field have been successful only for a dozen 
of absorbers at z ^ 1 (Le Brun et al. 1997, Turnshek et 
al. 2001; see also Kanekar & Chengalur 2003), but not 
in the redshift interval where DLAs are most commonly 
identified and investigated (z ^ 2). 

So far, the only method to estimate the dust abun- 
dance of individual DLAs relies on the comparison of the 
elemental abundances of refractory and non-refractory el- 
ements. Since the pioneer investigation on Cr and Zn by 
Pettini et al. (1994), the evidence for differential deple- 
tion of refractory elements relative to the volatile element 
Zn has accumulated over the years and is now well es- 
tablished (Hou et al. 2001, Prochaska & Wolfe 2002). An 
additional argument supporting the connection between 
depletion and dust is the correlation between depletion 
and fraction of molecular hydrogen (Levshakov et al. 2000, 
Ledoux et al. 2003). 



The existence of a correlation between differential de- 
pletion and metallicity, reported by Ledoux et al. (2003), 
suggests that the dust content of DLA systems evolves 
in the course of chemical evolution. However, in order to 
study the evolutionary properties of the dust we need to 
convert the differential depletions into dust-to-gas or dust- 
to-metal ratios. In a previous work we presented a method 
for estimating these ratios in DLA systems, taking also 
into account the depletion of Zn self-consistently (Vladilo 
1998; Paper I). In that work we derived a Z — V corre- 
lation between the metallicity corrected for dust effects 
and the dust-to-gas ratio. These two quantities, however, 
are not independent when the observed metallicity is very 
low. A much safer approach consists in studying the evo- 
lution of the dust-to-metal ratio, rather than that of the 
dust-to-gas ratio. The former is derived from the column 
densities of two metals and is less affected by error propa- 
gation than the latter, which requires, in addition, the use 
of the H I column density. 



A systematic study of the evolution of the dust-to- 
metal ratio in DLA systems is still missing. This type of 
study requires converting the differential depletions mea- 
sured at face value into the fraction of metals in dust 
form, fd, taking into account (i) variations of the dust 
composition that may occur in different types of galac- 
tic environments (e.g. galaxies with different metallici- 
ties) and (ii) nucleosynthetic effects that may affect the 
volatile/refratory abundance ratio (e.g. evolution of the 
Zn/Fe ratio). In a recent work we presented a method- 
ology aimed at tackling these problems (Vladilo 2002a; 
Paper II). This methodology was successfully applied to 
recover (i) the (stellar) abundance pattern of the Small 
Magellanic Cloud starting from interstellar measurements 
affected by dust depletion; (ii) the a/Fe ratio in DLA sys- 
tems corrected for dust effects (Vladilo 2002b; Paper III). 



In the present paper we apply this new method to de- 
rive the fraction of iron in dust form, /p e , for the sample 
of DLA systems with measurements of the Zn/Fe ratio. 
This sample now includes 38 absorbers, more than twice 
the sample of Paper I. We use the results of the most 
recent studies of Zn abundance in metal-poor stars to de- 
rive an educated guess for the intrinsic Zn/Fe ratio in DLA 
systems which, at variance with Paper I, is now a free pa- 
rameter of the method. Finally, for the first time in this 
type of investigation, the robustness of the results is tested 
against possible bias introduced by the observations or by 
the method itself. The paper is organized as follows. In 
the next Section we present the method of analysis and 
the sample. In Section 3 we present the results and test 
their robustness. In Section 4 we discuss the implications 
of the results in terms of formation and evolution of dust 
in primeval galaxies. The conclusions are summarized in 
Section 5. 
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2. Analysis of the observational data 

2.1. The method 

In order to study the evolutionary properties of dust in 
DLA systems we measure the fraction of metals in dust 
form in each system. We call dust fraction of an element 
X the number ratio 



/x = iV dust (X)/iV tot (X) , 



(1) 



where iVd us t (X) and N to t (X) are the column densities of X 
(atoms cm~ 2 ) in dust form and in the medium (gas plus 
dust), respectively. The dust fraction of a metal indicates 
how efficient the ISM of a galaxy is to incorporate part 
of its heavy elements into dust grains and represents, in 
practice, a dust-to-metal ratio by number. Here we adopt 
the dust fraction of iron, /p e , as an estimator of the dust- 
to-metals ratio in DLA systems. Iron is a natural choice 
since it is easy to measure in DLAs and is a reference 
element used in different types of abundance studies. Most 
important, iron is a refractory element which traces the 
dust even in the harshest interstellar conditions, where 
dust grains tend to be destroyed. This is demonstrated 
by the fact that in the warm halo gas of the Milky Way, 
where the level of depletions is known to be generally low, 
the dust fraction of iron is still quite high, fFe,wh — 0.77 
(Savage & Sembach 1996). 

2.1.1. Derivation of the dust fraction of iron 

The dust fraction of iron can be derived by comparing the 
observed column density of iron with that of a volatile el- 
ement X if we can make an educated guess of the intrinsic 
abundance ratio X/Fe in the DLA system. The method 
works best when the differential dust depletion between 
Fe and X is high. The most natural choice for X is Zn, a 
volatile element that tracks Fe well in the course of nucle- 
osynthetic evolution (see Section 2.2.3). If we call [Zn/Fe] s 
the educated guess of the intrinsic ratio in DLAs, we can 
use the equation 



/zn,,:10 [ i 



,(l+»7z„) + 
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101SN 



= o 
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derived in Paper III to obtain the normalized dust-to- 
metal ratio g = /fo//fc,4 from the ratio [Zn/Fe] measured 
in the DLA. 1 In the above equation the fractions in dust 
f-p e ,i and /zn,i are measured in one particular phase i of 
the Galactic ISM (e.g. cold gas, or warm gas in the disk or 
in the halo) . Without loss of generality, we adopt here the 
warm gas in the disk as a reference value: we stress that, 
at variance with Paper I, the present method allows for 
variations of the dust composition among different inter- 
stellar phases (see Paper II). These variations are taken 

1 [Zn/Fe] ee log{Af(Zn)/iV(Fe)}-log(Zn/Fe)Q, where iV(Zn) 
and iV(Fe) are the column densities derived from spectro- 
scopic measurements. The reference solar abundance by num- 
ber, log(Zn/Fe)Q = —2.86, is taken from meteoritic measure- 
ments (Anders & Grevesse 1989). 



into account by the parameter r]z n , which describes how 
the Zn/Fe ratio in the dust varies as a function of g. The 
parameter rj^ n is determined from studies of the Galactic 
depletion patterns (Paper II) . The method also allows the 
dust composition to vary as a function of the composition 
of the medium. This dependence is accounted for by the 
parameter ezn, discussed below. Finally, we stress that Eq. 
Intakes into account in a self-consistent way the fact that 
also Zn is partly depleted, while dust depletion studies in 
literature ignore this effect. 

Eq. FJ] is a transcendent equation in the unknown g 
that we solve by iterations. The equation has a unique 
solution inside the interval < g < Qu m , where gu m is the 
maximum possible value of g, corresponding to a fraction 
in dust of iron of 100%. Solutions outside this interval are 
not physical. 

2.1.2. Derivation of the metallicity 

Once the fraction in dust of iron is known, it is straightfor- 
ward to derive the total column density of iron (gas plus 
dust) and from this the iron abundance corrected for dust 
effects, [Fe/H] c , used here as the indicator of metallicity. In 
most cases the dust-corrected metallicity [Fe/H] c is very 
similar to the zinc abundance [Zn/H] taken at face value, 
which is commonly adopted in literature as the indicator 
of metallicity in DLA systems, with a typical difference 
[Fe/H] c - [Zn/H] < +0.1 dex. 

2.2. The parameters 

2.2.1. The dust fraction of zinc in the Galactic ISM 

Local interstellar studies indicate that a fraction of zinc 
is incorporated in dust form even in lines of sight with 
moderate level of depletion. In Paper II we derived the 
fraction in dust of zinc in the Galactic warm disk gas, 
fzn,wd — 0.59, using the same lines of sight with accurate 
HST measurements adopted by Savage & Sembach (1996) 
in their study of depletion patterns. Unfortunately, that 
sample of lines of sight included only a few zinc determi- 
nations. For this reason, we now derive fzn,wd also using 
the largest compilation of HST measurements of zinc and 
chromium (Roth & Blades 1995). As a result, we obtain 
fzn,wd = 0.30 using the same scaling law of interstellar de- 
pletions adopted in Paper II and imposing that the dust- 
corrected abundance of Cr in the local ISM is solar. We 
consider that the range 0.30 < fz n ,wd < 0.59 is fairly rep- 
resentative of the uncertainty of this parameter. 



2.2.2. The parameter ezn 

The parameter ezn describes how the Zn/Fe ratio in the 
dust composition varies as a function of the same ratio in 
the composition of the medium. By definition, ezn = 1 if a 
percent variation of the Zn/Fe ratio in the medium yields 
the same percent variation in the dust (Paper II). If, on 
the other hand, the Zn/Fe ratio in the dust is totally in- 
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dependent of variations of the ratio in the medium, then 
£zn = 0. In principle, £zn can be derived studying deple- 
tion patterns in galaxies with a known chemical composi- 
tion. A preliminary study of SMC data is consistent with 
£zn — 1- Here we consider both cases £z n = 1 and £zn = 
as extreme possibilities to solve Eq. [21 starting from ref- 
erence Galactic ISM values. It is important to note that, 
given the form of Eq. [21 in which £z n multiplies [Zn/Fe] g , 
the choice of ezn is not critical if the intrinsic Zn/Fe ratio 
is approximately solar. Luckily, the Zn/Fe is indeed close 
to solar and this explains why the results presented below 
are rather independent of the choice of £z n - 

2.2.3. The Zn/Fe ratio in DLA systems 

In order to make an educated guess of the Zn/Fe ra- 
tio in DLA systems we must rely on abundance studies 
of metal-poor stars. The early studies of Galactic stars 
yielded a Zn/Fe ratio close to solar (Sneden et al. 1991), 
but some recent work indicate that the Zn/Fe ratio can be 
overabundant at very low metallicity (Primas et al. 2000; 
Prochaska & Wolfe, 2002). The latter results have been 
used to support the notion that Zn/Fe ratio may decrease 
in the course of chemical evolution owing to a different 
nucleosynthetic origin of Zn and Fe (see e.g. Umeda & 
Nomoto, 2002). 

A large amount of Zn abundance measurements have 
been recently published by three different groups, all 
indicating that the Zn/Fe ratio is very close to solar, 
or slightly enhanced, in the range of metallicities typi- 
cal of DLA systems (Mishenina et al. 2002, Gratton et 
al. 2003, Nissen et al. 2003). The 87 measurements per- 
formed by Mishenina et al. , derived from spectra collected 
at the 1.93-m telescope of Haute Provence, yield a mean 
value <[Zn/Fe]> = +0.08 ± 0.13 dex (la) in the interval 
—3 < [Fc/H] < —0.2. The 48 measurements obtained from 
UVES/VLT spectra by Gratton et al. yield <[Zn/Fe]> 
= +0.13 ± 0.13 dex in the interval -2.4 < [Fe/H] < -0.3. 
The 29 measurements by Nissen et al., also obtained with 
UVES/VLT data, yield <[Zn/Fe]> = +0.04 ± 0.08 dex in 
the interval -2.4 < [Fe/H] < -0.7. The scales of the zinc 
abundances in these measurements might be in error up to 
0.1 dex owing to differences between the photospheric and 
meteoritic solar abundance. For instance, Gratton et al. 
used a reference solar ratio log(Zn/Fe)Q = —2.95 derived 
from a consistent analysis of the solar spectrum. Using 
instead the solar meteoritic value log(Zn/Fe)Q = —2.86 
(Anders & Grevesse 1989), one would derive <[Zn/Fe]> 
= +0.04 ± 0.13 dex from Gratton et al.'s data. Also a 
systematic difference in the scale of effective temperature 
could easily account for the mean offset relative to the 
solar value. 

The [Zn/Fe] data do not show a trend with metallic- 
ity, with the possible exception of a very weak decrease 
with [Fe/H] in the Nissen et al.'s sample and in the sub- 
sample of halo stars of Mishenina et al., neither statisti- 
cally significant. Therefore a constant value of [Zn/Fe] g 



is a reasonable choice in the metallicity interval typi- 
cal of DLA systems. We considered two possible values, 
namely [Zn/Fe] g = 0, representative of the solar value, 
and [Zn/Fc]g = +0.1, representative of a modest enhance- 
ment. The potential effects of an hypothetic decrease of 
Zn/Fe with Fe/H are considered in the discussion. 



2.3. The sample 

The data used in the present investigation were selected 
searching for DLA systems with published measurements 
of Fc II and Zn II column densities. All systems with N(H i) 
in excess of the canonical threshold 7V(Hi) = 10 20 3 atoms 
cm~ 2 , or consistent with the threshold at ~ lcr level, were 
included in the sample. Upper limits were excluded. The 
resulting sample, listed in Table 1, includes 38 systems 
in the redshift interval 0.6 < z < 3.4. The references to 
the original measurements are given in the table. Most 
of the measurements were obtained using high resolution 
spectrographs fed by 8-10m class telescopes such as the 
Keck and the VLT. Care has been taken in checking that 
all the column densities are derived with an updated and 
consistent set of oscillator strengths (Bergeson & Lawler 
1993; Welty et al. 1999). 

3. Results 

In the last 4 columns of Table 1 we list the fraction in 
dust of iron in DLA systems /f c = Q f¥c,wd obtained 
from Eq. [21 for different choices of the input parameters 
fzn,wd, £zn and [Zn/Fe]g, specified in the notes to the ta- 
ble. The quoted errors have been derived by propagat- 
ing the Fen and Znii column-density errors. The results 
obtained for a solar [Zn/Fe] g are labeled S; those for an 
enhanced [Zn/Fe] g (= +0.1) are labeled E. For the case 
S we list the results obtained using fz a ,wd = 0.59 and 
fzn,wd = 0.30, labeled Sa and Sb, respectively. The differ- 
ences are in most cases well below the quoted errors. The 
results S do not depend on the adopted value of £z n (see 
Eq. [21 . The results E do depend on this parameter and 
are labeled E0 and El to indicate the cases £z n = and 
= 1, respectively. Comparing the columns E0 and El of 
Table 1 one can see that the dependence on £z n is indeed 
weak, as anticipated above. 

3.1. Frequency distribution of depletions 

As a first step of our analysis we consider the frequency 
distribution of iron depletions 2 , 

5 Fc = log(l - / Fe ) , (3) 

2 It is easy to show that Eq. [H] is equivalent to the clas- 
sical definition S Fc = log[V(Fe) b s /V(H) ob s] - log(Fe/H) , 
when the intrinsic abundance of the medium is solar, as in 
the Galactic ISM. 
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a The fraction in dust of iron is derived from the expression I2I and represents a dust-to-mctal ratio by number. Ii 
QSO 0454+039 the measured [Zn/Fc] ratio is lower than the adopted [Zn/Fc] g for the parameter sets E0 and El; 
avoid a negative, unphysical result. For the same reason, the error bars of /fg have been truncated in a few cases. 
b Logarithm of the observed column densities expressed in atoms cm -2 . 

c Parameter set Sa: [Zn/Fc] g — 0.0 dex, fzn,wd — 0.59; results independent of the adopted value of £zn- 
[Zn/Fe]g — 0.0 dex, fzn,wd — 0.30; results independent of the adopted value of ezn- 
[Zn/Fe] s = +0.1 dex, / Zn-rod = 0.59; e Zn = 0. 
[Zn/Fc] K = +0.1 dex, f Z n, w d = 0.59; e Zn = 1. 



1 the z — 0.8597 system towards 
in these cases we set /f c — to 



Parameter set Sb 
Parameter set E0 
Parameter set El 



Sum of all components of the Fc II and Z11 II column densities. 
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shown in Fig. ^ For the sake of comparison with Galactic 
interstellar studies we also plot the depletions represen- 
tative of the Galactic cold disk gas (dashed-dotted line), 
the warm disk gas (dashed line), and the warm halo gas 
(dotted line), taken from Savage & Sembach (1996). 

The frequency distributions are shown for different 
choices of input parameters (cases Sa, Sb and El). The 
main results common to all cases can be summarized as 
follows. The relative number of systems decreases with in- 
creasing depletion |<5f |- A large fraction of systems have 
depletions lower than the lowest Galactic depletions (those 



of the warm halo gas; dotted line). The remaining DLAs 
have depletions similar to those of the warm halo or disk 
gas; only for these cases it is correct to conclude that 
"depletions in DLA systems are typical of warm Galactic 
gas , a claim often quoted in literature. Cold disk deple- 
tions are not found in the sample. These conclusions are 
also valid for the case EO, not shown in the figure. 
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Fig. 1. Distribution of iron depletions for the DLA systems in Table 1. Panels a), b) and c) correspond to cases Sa, Sb 
and El, respectively. Vertical lines: typical values of interstellar iron depletions in the Milky Way (dotted line: warm 
halo gas; dashed line: warm disk gas; dotted-dashed line: cold disk gas). 




[Fe/H]„ 




Fe/H] c 




Fig. 2. Iron dust fraction versus metallicity for the sample of DLA systems listed in Table 1. Panels a), b) and c) 
correspond to cases Sa, Sb and El, respectively. Horizontal lines: representative values of iron dust fractions in the 
Galactic ISM (legenda as in Fig.^J. Solid curve: empirical law (@J describing the evolution of the dust fraction with 
metallicity. 



3.2. Dust fraction versus metallicity 

We next consider the relation between dust fraction, f-p e , 
and metallicity, which is shown in Fig.[21for three different 
choices of input parameters (cases Sa, Sb and El). A very 
similar relation is found for the case EO, not shown in the 
figure. The main result, common to all cases considered, 
is that the dust fraction increases with metallicity, with a 
flattening at [Fe/H] c > — 1 dex and a significant scatter. 
The data follow an empirical law of the type 

he « /Fc,o ~ 10«« Fc / H MFe/HU (4) 

with / Pe>0 ~ 0.9, a ~ -1.2, and [Fe/H] D ~ -1.7 (solid 
curve in the figures). The upper (lower) envelope of the 
data, shown in Fig. |3Jd, can be represented with the same 
law, with parameters /Fe,o — 0.99 (0.67), a ~ — 1 (—4), 
and [Fe/H] D ~ -1.7 (-2.15). 

The comparison with local ISM depletions reveals that 
the systems with depletions lower than in the Galactic 
halo gas are mostly found at [Fc/H] c < — 1 dex. 



We stress that the increase of /p e with metallicity is 
confirmed even assuming that [Zn/Fe] g evolves from m 
+0.2/ + 0.3 dex at [Fe/H] < - 2 dex, down to ~ at 
solar metallicity, as claimed by some authors (see Section 
2.2.3). In fact, in this case we would obtain even less dust 
than in Fig. 2c at [Fe/H] < — 2 dex, but the same amount 
of dust as in Fig. 2a at solar metallicity. One can see that, 
as a consequence, the rise of /Fe with metallicity would be 
even steeper. We now discuss the stability of the result in 
light of possible observational bias. 

3.3. Selection bias 

Several types of observational bias are known to affect 
the measurements of the AT(Hi), iV(Feii), and iV(Znil) 
column densities. For instance, the detection limit 3 of 

3 The detection limits are estimated at 3a level for a signal- 
to-noise ratio S/N ~ 30 and resolving power R = A/(AA) ~ 
5 x 10 4 at A ~ 500 nm. Oscillator strengths for the Znll 
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the strongest line of the Zn n resonance doublet prevents 
the measurement of column densities A(Znii) ^ 10 11 ' 75 
atoms cm -2 . The detection limit of Fen depends on 
the wavelength coverage of the spectra, since a variety 
of strong Fen transitions are present at different wave- 
lengths. Assuming that we do not cover the stronger tran- 
sitions at longer rest wavelengths, we obtain a conserva- 
tive limit from the detection limit of the 160.8 nm line, 
A(Feii) xo 12 - 88 atoms cm" 2 . Also Hi measurements are 
affected by selection bias. The range of 7V(H i) is limited 
on one side by the definition threshold in DLA systems 
A(H i) > 10 20 3 atoms cm -2 and, on the other side, by the 
fact that values above iV(Hi) > 10 20 - 85 atoms cm" 2 have 
never been observed, even if they are detectable. Finally, 
the obscuration of the background QSO by dust extinc- 
tion in the DLA system could be responsible for the lack 
of systems with high metal column densities, namely those 
with A(Znn) ^ 10 13 - 1 atoms cm" 2 (Boisse et al. 1998), 
which are also detectable. 

In order to assess the effect of these observational bias 
on the present study we built a mock sample of n DLA 
systems with a flat distribution of metallicities and de- 
pletions (or dust fractions) . We then assigned to each sys- 
tem a random-generated H I column density in the allowed 
range 20.3 < logA(Hi) < 20.85 and computed A(Feii) 
and TV(Znil) for each pair of ([Fe/H],/p ) values. Finally, 
we rejected cases with A(Znii) and TV(Feii) outside the 
limits discussed above. At this point, we applied to the 
biased sample the same method for the derivation of dust 
fractions and metallicities described in Section 2.1 and 
the same analysis presented in Sections 3.1 and 3.2. In 
this way, we tested the potential effects of the column- 
density bias and, at the same time, of the application 
of the method itself. The results are shown in Fig. for 
n = 1000 systems (before rejection). 

The frequency distribution of depletions of the biased 
mock sample is shown in the panel (a) of the figure. One 
can see that, in spite of the large fraction of data points 
rejected (65%), the distribution is still approximately flat 
in the range —2.2 < S(Fc) < 0, as defined by construc- 
tion. This indicates that the overall distribution of deple- 
tions is not significantly affected by bias and, therefore, 
the steepness of the real distribution in Fig. ^ is genuine. 
We conclude that the deficiency of systems with high levels 
of dust depletion seems to be an intrinsic property of DLA 
systems rather than a consequence of selection effects. 

The plot of dust fraction versus metallicity, for the 
mock sample built to uniformly cover the range < fp c < 
1, is shown in Fig. |3Jd. The non-rejected data of the mock 
sample (crosses in the figure) are able to populate regions 
outside the boundaries in which the real data are con- 
fined (solid curves). This indicates that the selection bias 
is not responsible for the lack of real data outside these 
boundaries. In addition, the rather homogeneous distribu- 

202.6136 nm and the Fell 160.8451 nm lines are taken from 
Bergeson & Lawler (1993) and from Welty et al. (1999), re- 
spectively. 
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Fig. 3. Test of selection effects performed on a mock sam- 
ple of n = 1000 DLA systems with Hi, Fen and Znii 
column densities rejected as explained in Section 3.3. a) 
Frequency distribution of iron depletions resulting (after 
rejection) for an adopted flat distribution of <5(Fe). b) Dust 
fraction versus metallicity resulting for an adopted flat dis- 
tribution of fp e ; the solid curves represent the envelopes of 
the real data in Fig.0 The depletions of the mock sample 
have been modeled with the set of parameters Sa (Table 
1); similar results are obtained for the other cases. 

tion of the data in Fig.^, derived from the application of 
our method to the mock column densities, indicates that 
the method itself does not induce artificial trends between 
/Fe and [Fe/H]. We conclude that the rise of dust fraction 
with increasing metallicity is a genuine property of DLA 
systems. 

Clearly, a better understanding of the dust obscura- 
tion effect is required in order to take into account the 
role of selection bias in a more realistic way. In any case, 
we expect the extinction of the QSO to be more and more 
effective with increasing metallicity of the intervening ab- 
sorbers. Therefore, the fact that DLA systems with high 
dust content (/Fe > 0.8) are detected at high metallicity 
demonstrates that the obscuration bias cannot be respon- 
sible for the lack of such systems at low metallicity (com- 
pare the data points with [Fe/H] < — 1.5dexand > —1.5 
dex in Fig.|5J). On the other hand, the number of absorbers 
with depletions typical of cold interstellar gas and high 
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metallicity can be affected by dust obscuration. Therefore 
the deficit of cold gas depletions could be less severe than 
shown in Fig. ^ but only for the systems of high metal- 
licity. 

3.4. Time evolution 

Finally we study the time evolution of the dust fraction 
after converting the redshift of each absorber to the look- 
back time or age of the Universe 4 . The advantage of this 
approach is that we compare /f c with a completely inde- 
pendent observable, while in the study versus metallicity 
we compare quantities derived from the measurements of 
the same column densities. However, the interpretation of 
the time evolution is complicated by the fact that DLA 
galaxies may start their chemical evolution at different 
cosmic times. Inspection of the plot of /p c versus age of 
the Universe, shown in Fig. 0^,, does not reveal a clear 
trend, even though the data are consistent with a sudden 
rise of the dust-to-metal ratio between w 2 and 4 Gyr af- 
ter the Big Bang. Some evidence of evolution is found by 
binning all the systems seen before and after the median 
value of look-back time of the present sample, t\ ~ 10.5 
Gyr, and comparing the frequency distribution of their de- 
pletions (Fig. HJ>)- O ne can see that the two distributions 
are slightly shifted, the DLAs seen at the earliest age of 
the universe (shaded histogram) having lower depletions 
|<5(Fe)| than those observed in more recent times (empty 
hystogram) . Again, these results hold valid for all the pos- 
sible choices of input parameters (case Sa, shown in the 
figure, and cases Sb, E0 and E, not shown). 



4. Discussion 

We now discuss the empirical results presented in the pre- 
vious section assuming that they are substantially un- 
affected by selection bias. We first try to understand 
whether the observed trends are mainly governed by vari- 
ations of the physical properties or by chemical evolution. 
An example of the difficulty of disentangling these two ef- 
fects is given in Fig. [IJ,, where we plot the iron depletions 
versus l c oc 7V(C u*)/N(H i), an indicator of the cooling 
rate of the gas (Wolfe et al. 2003a). Apparently, |<5(Fe)| 
seems to increase with log l c and this trend could be taken 
as evidence for a correlation of the depletions with some 
physical property of the gas. However, l c is by definition 
proportional to the metallicity and the trend in Fig. [S^, 
may simply mirror the rise of /F e with metallicity shown 
in Fig. We now consider other observational evidence 
that might help in disentangling the role of physical pro- 
cesses. 



4 We calculated these conversions for a flat universe with 
H = 71 km s" 1 Mpc" 1 and n m = 0.27 (Spergel et al. 2003). 



4.1. The dependence on physical parameters 

In the Galactic ISM, each interstellar phase is character- 
ized by a particular level of depletion. Cold, high-density 
gas, where the dust has higher accretion probability and 
is more shielded by destruction processes, is character- 
ized by high depletions; warm, low-density gas is instead 
known to have low depletions (e.g. Savage & Sembach 
1996). This link between the physical state and level of 
depletion might be at work also in the ISM of other galax- 
ies, including DLA systems. If so, low values of depletions 
would map regions of warm gas, while higher depletions 
colder regions. According to this interpretation, the fre- 
quency distribution in Fig. ^ would indicate that in DLA 
systems warm regions are more frequent (have a higher 
cross-section) than cold regions. The evolutionary trends 
with metallicity (Fig. EJ) and cosmic time (Fig. 0J would 
indicate that cold regions are extremely rare in the early 
stages of evolution, but eventually they become more fre- 
quent. 

Some evidence in favour of this interpretation comes 
from measurements of the spin temperature, T s (K), de- 
rived 21-cm absorption observations (Carilli et al. 1996, 
Chengalur & Kanekar 2000, Kanekar & Chengalur 2003) 
and from studies of the AT(Cll)/JV(Cl) ratio (Liszt 2002), 
which suggest that the warm phase is generally dominant 
in DLA systems, especially at high redshift (see however 
Wolfe et al. 2003b). Nevertheless, a detailed study of the 
results on a system-by-system basis does not support a 
direct connection between depletions and physical param- 
eters in DLAs. In fact, if the observed trends of depletions 
were governed by variations of the physical properties, we 
would expect to find a decrease of |<5(Fe)| with increas- 
ing T s or A r (Cn)/7V(Ci), i.e. with increasing importance 
of the warm component. The comparison of depletions 
versus T s for the systems in common with the compila- 
tion by Kanekar & Chengalur (2003) does not support 
this expectation (Fig. Et>). Also the comparison with the 
N(C li)/N(C i) determinations by Liszt (2002) does not 
reveal the expected trend (Fig. Et)- 

We also compared the depletions with the molecular 
fraction /(H 2 ) = 2Ar(H 2 )/[Ar(Hi)-t-2Ar(H 2 )], which is sen- 
sitive to physical conditions. Some evidence for a trend be- 
tween differential depletions and /(H 2 ) has been reported 
in previous work (Levshakov et al. 2000; Ledoux et al. 
2003). The limited number of measurements in common 
with the present sample are marginally consistent with the 
existence of such a trend (Fig. Oi) . In any case, this trend 
would only confirm the key role of dust in the process of 
H 2 formation, but would not provide evidence for a direct 
link between the physical properties and the abundance of 
dust. In fact, the correlation is expected because dust acts 
as a catalyst of H 2 formation and provides a shield against 
photo-dissociation of molecules by ultraviolet radiation. 

In summary, there is no evidence that the evolution 
of the dust fractions/depletions is primarily governed by 
variations of the physical state of the gas. In the follow- 
ing we consider the chemical evolution as the main driver 
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Age of the Universe (Gyr) 5(Fe) 

Fig. 4. Time evolution of the dust content in DLA systems, a) Diamonds: /p e versus age of the universe; dashed area: 
reionization epoch 6 < z < 20 (Fan et al. 2003; Spergel et al. 2003); solid, dotted and dashed curves: Inoue's (2003) 
model predictions of the evolution of the relative dust-to-metal ratio for different sets of input parameters (see text 
for explanations); the curves have been normalized to the present-day value and scaled assuming a formation redshift 
Zf = 20. b) Frequency distribution of iron depletions in two bins of look-back time (shaded area: ti > 10.5 Gyr; empty 
histogram: ti < 10.5 Gyr). The fraction in dust are computed for the set of parameters Sa. 



of the evolution of the dust abundance in DLA systems. 
However, we expect that variations of the physical state 
of the gas will contribute to the scatter of dust abundance 
observed at a given metallicity. 

4.2. The dependence on chemical evolution 

The regular trend shown by the data in Fig. [3 is an ar- 
gument in favour of a direct dependence of depletions on 
the metallicity, i.e. on the chemical evolution. 

The rise of the dust-to-metal ratio fp e in the course of 
chemical evolution is unexpected, given the approximately 
constant dust-to- metal ratio V/Z measured in local galax- 
ies (Issa et al. 1990; Schmidt & Boiler 1993). However, the 
two results are not in contradiction since the data in Fig. 
suggest the existence of a flattening of fp c rs V/Z at 
[Fe/H] > — 1 dcx, i.e. when the metallicity approaches 
values typical of the local sample of T>/Z measurements. 

Also the scatter of /f c in DLA systems is consistent 
with the observed properties of the V/Z ratio in present- 
day galaxies. In fact, a spread of the T>/Z ratio is found 
when different types of galaxies are compared, such as the 
blue compact dwarfs and dwarf irregulars investigated by 
Lisenfeld & Ferrara (1998). Therefore, the scatter of fp e 
may simply reflect the inhomogeneity of the population of 
DLA systems, which in fact is believed to include galaxies 
of different types. The V/Z ratio about 1/4 to 1/7 of 
the Galactic one recently found in a dwarf galaxy with 
an average metallicity ~ 1/4 solar (Galliano et al. 2003) 
is a further example indicating that the properties of the 
dust fraction fp c in DLA systems are consistent with the 
properties of the dust-to- metal ratios V/Z in the local 
universe. Thefore, the trend seen in Fig. [5] may represent 
a general characteristic of galactic evolution, rather than 
a specific property of DLA galaxies. 



The rise of V/Z with metallicity suggests that the ef- 
ficiency of dust formation is very low at the earliest stages 
of chemical evolution but increases regularly as metals be- 
come more abundant. As a first step in interpreting this 
effect, we consider the possibility that the most common 
processes of dust formation may show a dependence on 
the metallicity. Dust can be formed in the ejecta of type 
II SNe (Moseley et al. 1989; Elmhamdi et al. 2003; Dunne 
et al. 2003) and in the cool winds of late-type giant stars 
(see Sedlmayr 1989 and refs. therein). It can probably also 
be formed in the ejecta of type la SNe, even though in 
this case the evidence is not direct (Clayton et al. 1997). 
Other processes of dust formation are believed to give 
a minor contribution on a galactic scale (see e.g. Dwek 
1998). The composition of the dust is poorly known on 
observational grounds and often inferred from chemical 
equilibrium condensation calculations. The relative abun- 
dance of iron in dust grains is not known, even though 
iron could be present both in carbon-rich dust (e.g. FesC; 
Lattimer et al. 1978), and in silicate grains (e.g. FeSiOa, 
Fe 2 Si0 4 ; Ossenkopf et al. 1992); the dust formed by SNe 
la, if any, is expected to be iron-rich. In the following we 
consider the possible dependence on the metallicity of the 
main mechanisms of dust formation, assuming that a sig- 
nificant fraction of iron can be incorporated in dust form 
in such processes. 

Given the short life-times of their massive progenitors, 
SNe II probably represent the only source of dust at very 
high redshift. In fact, the time-scale for dust production 
by the winds of red giants or by SNe la is of the order of 1 
Gyr, i.e. larger than the age of the universe at z 4. The 
existence of dust at z > 4 (e.g. Carilli et al. 2000) is a strin- 
gent argument in favour of the production of dust by SNe 
II in the early universe. The systems with highest redshift 
in our sample are seen at about JS; 1.8 Gyr after the start 
of the reionization epoch (Fig. 0]i) and a significant part 
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of their dust must have been produced by SNe II. A dust 
formation mechanism for SNe II dependent on metallicity 
has been found by Todini & Ferrara (2001). In this pro- 
cess the density of heavy elements in the ejecta becomes 
large enough at higher metallicities to allow the state of 
supersaturation to be reached more easily, favouring the 
production of dust. The metallicity-dependent dust forma- 
tion in SNe II found by Todini & Ferrara, together with 
the lack of other sources of dust, may explain the very low 
dust fraction typical of most of the highest-redshift sys- 
tems in our sample (shaded histogram in Fig.^J}), which 
have the lowest metallicities. 

After w 1 Gyr from the start of the chemical evolu- 
tion the dust starts to be injected also by the winds of 
red giants and possibly by SNe la. We consider three pos- 
sibilities: (1) these additional contributions are negligible 
compared to the dust production of SNe II; (2) the pro- 
duction of dust in the wind of late- type giants is dominant; 
and (3) the production by SN la is dominant. 

In the first case the metallicity-dependent dust pro- 
duction of SNe II may explain the rise of the dust fraction 
with metallicity not only at the highest rcdshifts, but also 
at later stages of evolution. 



In the second case another mechanism may provide a 
metallicity-dependent dust formation, namely the depen- 
dence on the metallicity of the wind efficiency. In fact, re- 
cent work on dust production in giant stars indicates that 
the efficiency of the stellar winds is lower at lower metal- 
licity (Ferrarotti & Gail 2003). Therefore, the gradual rise 
of the efficiency of their winds with increasing metallicity 
may yield a gradual rise of the dust fraction. However, 
this mechanism can explain the rise of /F e only if iron is 
already present in the atmosphere of the late-type giants, 
since iron is not synthesized by these stars. This implies 
that this mechanism could only work at later stages of 
evolution, ps 1 Gyr after iron is already present in the gas 
out of which the future late-giants are born. 

In the third case, the (hypothetical) iron-rich dust from 
SN la could provide an additional interpretation: the time 
delay between the early production of dust by SNe II and 
the subsequent production of iron-rich dust by SNe la 
would work as a mechanism for increasing the iron dust 
fraction with time. A potential problem with this explana- 
tion is the possibility that SNe la can only explode when 
[Fe/H] ^ -1 dex (Kobayashi et al. 1998), i.e. after the 
rapid increase of /F e seen in Fig [2 



G. Vladilo: The early build-up of dust in galaxies 



11 



In summary, all the important mechanisms of dust for- 
mation might be able to provide a rise of the dust fraction 
in the course of chemical evolution. The most convinc- 
ing mechanism, however, is the metallicity-dependent pro- 
duction by SNe II proposed by Todini & Ferrara (2001). 
Clearly, the evolution of the dust content will be deter- 
mined not only by the dust formation rate, but also by 
the destruction and accretion rates, not discussed here 
in detail. By assuming that the ratio of the destruction 
and accretion time-scales is proportional to the time-scale 
of star formation, Inoue (2003) predicts that the dust- 
to-metal ratio should increase in the course of chemical 
evolution. The rise of the dust-to-metal ratio predicted 
by the "standard Galactic model" of Inoue (solid curve 
in Fig. 0k,) is less steep than the rise observed in DLA 
systems. The better agreement with the case of "no SNe 
destruction" (dotted line) may simply reflect the fact that 
iron is one of the most refractory elements, not easily de- 
stroyed by the SNe shocks. The range of dust-to-metal 
ratios predicted by the three models shown in the Fig. 6 
of Inoue (2003) approximately bracket the observed data. 
The highest values of /F e observed at ss 2 Gyr after the Big 
Bang suggest a very early epoch of galaxy formation for 
these particular DLA systems, since even the model with 
the highest accretion rate (dashed line) has some difficulty 
in reproducing these high values of /p e , even assuming a 
formation redshift zt = 20. 

The above conclusions must be considered tentative 
given the simplified nature of the analytical models and 
because the conversion between /p e and the dust-to- 
metal ratio by mass is model-dependent (e.g. the conver- 
sion requires some assumption on the dust composition). 
Nevertheless, this exercise is an example that shows how 
the present data can be used to constrain galactic evolu- 
tion models incorporating the cycle of dust formation and 
destruction. The implementation of the dust component 
in realistic models calibrated on galaxies of the local uni- 
verse can in turn be used to provide fresh clues on the 
redshift of formation and the nature of DLA galaxies (e.g. 
Calura et al. 2003). 

5. Conclusions 

We have investigated the evolutionary properties of the 
dust-to-metal ratio in DLA systems with the aim of cast- 
ing light on the early build-up of dust in galaxies. The 
results that we obtain are representative of high-redshift, 
metal-poor galaxies selected on the basis of their Hi ab- 
sorption column density, i.e. A(Hi) > 2 x 10 20 atoms 
cm~ 2 . 

As an estimator of the dust-to- metal ratio in DLA sys- 
tems we have used the fraction of iron atoms in dust form, 
/Fe- One advantage of this approach is that iron is a re- 
fractory element which traces the presence of dust even 
in harsh interstellar environments. The dust fraction /p e 
has been derived by comparing the iron abundance with 
the abundance of the volatile element Zn. As in Paper I, 
we take into account the fact that also Zn can be incor- 



porated in dust form. However, the methodology adopted 
here features important improvements over the one fol- 
lowed in Paper I: (1) the intrinsic Zn/Fe ratio in DLA 
systems is now a free input parameter; (2) the dust chem- 
ical composition is allowed to vary in different galactic 
environments (see Paper II); (3) we prefer the dust-to- 
metal ratio, rather than the less reliable dust-to-gas ratio, 
to study the evolution of the dust content. The main re- 
sults of the present study can be summarized as follows. 

The frequency distribution of iron depletions, <5(Fe) = 
log(l — /Fe), peaks at S(Fe) ~ 0, and decreases regularly 
with increasing |<5(Fe)|. A significant fraction of DLA sys- 
tems has depletions lower than the lowest values measured 
in the Galactic ISM, i.e. those of the Galactic halo gas. 
The frequently quoted claim that depletions of DLA sys- 
tems " are typical of Galactic warm gas" is correct only for 
the remaining systems. High values of depletion, typical 
of Galactic cold disk gas, are not found. 

The dust fraction /p c increases with metallicity, al- 
beit with a significant scatter and a flattening of the 
trend when [Fe/H] <J — 1 dex. An empirical relation be- 
tween dust fraction and metallicity of the type /f c ~ 
/ Fe , D - lO^/HMFVHU, with / Fc ,o 0.9, a ~ -1.2, 
and [Fe/H] ~ —1.7, provides a good fit to the data. 

The comparison between the dust fraction /p c and the 
look-back time, ti (derived from the absorption redshift), 
reveals a weak evidence for a rise of /f c with cosmic time. 
In spite of a large scatter, the data are consistent with a 
fast rise of /f c starting at about 2 Gyr after the Big Bang 
and lasting a few billion years. Indeed, the comparison of 
the frequency distributions of depletions for the two sub- 
samples with ti > 10.5 Gyr and t\ < 10.5 Gyr suggests 
that the dust fraction does increase with cosmic time. 

The above conclusions are rather stable for different 
choices of the input parameters required for the deter- 
mination of /f c . In particular, they are not affected by 
changes of the adopted guess of the intrinsic Zn/Fe ratio 
in DLA systems, [Zn/Fe] g , in line with the measurements 
of the same ratio in metal-poor stars. The large amount 
of measurements recently reported by three independent 
groups consistently indicate that the Zn/Fe ratio in metal- 
poor stars is very close to solar and probably constant in 
the metallicity interval typical of DLA systems. In any 
case, assuming that the intrinsic Zn/Fe ratio decreases 
with chemical evolution, as reported by some authors, the 
existence of the trends with metallicity and cosmic time 
would even be reinforced. 

We have tested the robustness of the results in light of 
the selection bias which are known to affect the measure- 
ments of the Hi, Fen, and Znii column densities. From 
the analysis of a mock sample of DLA systems with as- 
signed values of dust fraction and metallicity we conclude 
that the bias effects are small and therefore that the trend 
between dust fraction and metallicity represents a genuine 
property of DLA systems. This is also true for the dis- 
tribution of depletions, which is peaked at low values and 
shows a deficiency of systems with depletions typical of 
Galactic cold disk gas. The dust obscuration bias may in 
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part be responsible for such deficiency, but probably only 
for the systems of high metallicity. 

The presence of different interstellar phases with a 
range of physical parameters in DLA systems may con- 
tribute to the scatter of the dust fractions measured at 
each metallicity. However, the depletions do not show cor- 
relations with physical parameters measured in DLA sys- 
tems, such as the spin temperature or the Cn/Ci ioniza- 
tion ratio. We propose that the evolutionary trends of the 
dust fraction in DLA systems are mostly governed by the 
chemical evolution of the galaxies where the absorptions 
originate, rather than by the evolution of their physical 
properties. 

The extremely low dust fraction of iron measured at 
[Fe/H] < — 2 dex suggests that dust formation is very in- 
efficient during the earliest stages of galactic chemical evo- 
lution. We considered the possibility that the main mecha- 
nisms of dust formation in galaxies may have a metallicity- 
dependent efficiency. Indeed, previous work already indi- 
cates that the dust formation efficiency in the ejecta of 
SNc II and in the winds of late-type giants may increase 
with metallicity. We argue that, in addition, the injection 
of iron-rich dust by type SNe la may yield a rise of /f c 
in the course of evolution. The metallicity-dependent dust 
formation in SN II seems to be the most promising candi- 
date for explaining the rise of fee at [Fe/H] ^ — 1 dex. 

In order to reproduce the observed evolution of the 
dust content of DLA systems a simultaneous treatment 
of the main processes of dust formation, accretion and 
destruction is required. The observed increase of /f c with 
metallicity is contrary to the expectations of current chem- 
ical evolution models, which postulate or infer an approxi- 
mately constant dust-to-metal ratio, with the exception of 
the simple, analytical models by Inoue (2003), which pre- 
dict an increase of the dust-to-metal ratio D/Z, starting 
from very low values. The present results may represent 
a unique guideline for incorporating the dust component 
in well-tested models, particularly for the early stages of 
evolution. The successful implementation of the observed 
trends between /F e , metallicity and cosmic time in realistic 
models will provide fresh clues on the formation rcdshift 
and nature of DLA galaxies. 
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